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ABSTRACT 

By exploiting two sets of high-resolution images obtained with HST 
ACS/WFC over a baseline of ~ 10 years we have measured relative proper mo¬ 
tions of ~ 70, 000 stars in the stellar system Terzan 5. The results conhrm the 
membership of the three sub-populations with different iron abudances discovered 
in the system. The orbit of the system has been derived from a hrst estimate of its 
absolute proper motion, obtained by using bulge stars as reference. The results of 
the integration of this orbit within an axisymmetric Galactic model exclude any 
external accretion origin for this cluster. Terzan 5 is known to have chemistry 
similar to the Galactic bulge; our Endings support a kinematic link between the 
cluster and the bulge, further strengthening the possibility that Terzan 5 is the 
fossil remnant of one of the pristine clumps that originated the bulge. 

Subject headings: proper motions: general; stellar system: individual (Terzan 5); 


Based on observations (G012933, G09799) with the NASA/ESA Hubble Space Telescope, obtained at 
the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
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1. INTRODUCTION 


Terzan 5 is a stellar system located at the edge of th e inner bulge of th e Galaxy 


(/ = 3.8395°, b = 1.6868°) at a distance of 5.9 kpc from the Sun ( Valenti et ah 2007h. Becaus e 


of the large and spatially variable extinction in that region of the sky flMassari et al.ll2012l) . 
observations of this system are extremely challenging, particularly in the optical bands. This 
is why the complex natur e of Terzan 5 has b een revealed only rec ently by means of high 
resolution IR photometric (IFerraro et al.ll2009l) and spectroscopic fsee lOriglia et al.ll201ll) ob¬ 
servations of its stellar populations with the ESO-VLT and Keck telescopes. These studies 
revealed that Terzan 5 hosts two distinct populations with signihcantly different iron content 
(A[Fe/H]= 0.5 dex), different level of a-element enhancement, and no evidence of the anti- 
corre lation among light elements commonly observed in globular clusters (GCs; e.g.. lGarretta 


2014J and references therein). The radial distribution of the two populations is incompatible 


with that of held stars. In addition, the metal-rich compone nt has been found to be signih¬ 
cantly more centrally concentrated than the metal poor one flLanzoni et al.ll2010l). Also, the 


two p opulations have been found to share the same mean radial velocities (iMassari et ah 


2014al and references therein) and the same center of gravity flTvanzoni et al.ll201011 . All these 


features add weight to the membership of these populations. However, even with the spa¬ 
tial distribution and radial velocity memb ership secure, their actual membership has been 


questioned (e.g. IWillman fc Stradeiil2012l) because of the strong held star co ntamination af 


fecting many GGs in the direction of the bulge (see Valenti et ah 2007 . 201oh. Recently, als o 


a third, metal-poor and a-enhanced population has been discovered flOriglia et al.N2013l) . 


cluster (jOriglia e 

ah 

2013 

) and it has survived to all statistical decontamination tests (see 

the discussion in 

Massari et ah 

2ni4bh 


In order to properly separate genuine Terzan 5 stars from foreground and background 
sources, we measured high-precision relative proper motions (PMs) of individual stars in the 
direction of the sy stem, by exploiting the superb a s trometric accuracy of the Hubble Space 


Telescope (see 

e.g. 

Anderson & van der Marel 

2010; 

Bellini et ah 

2014; 

Watkins et ah 

2015; 

Massari et ah 

2013 

Dalessandro et ah 

2 OI 3 I. and references thereinh Here we report on the 


results of this proper motion study, which also conhrms the membership of the three iron 
groups. 


This paper is organized as follows. In Section |2] we present the used data set, while 
in Section |3] we describe the techniques adopted to measure PMs. Finally, we present the 
results of our work in Section H] and we summarize the conclusions in Section [5l 
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2. OBSERVATIONS AND DATA REDUCTION 


In the context of a multi-wavelength program aimed at studying the photometric proper¬ 
ties of the multi-iron sub-populations of Terzan 5, we have acquired (GO-12933, PI: Ferraro) 
a set of optical and near-infrared (IR) images of the system, by using the Wide Field Chan¬ 
nel (WFC) of the Advanced Camera for Surveys (ACS) and the IR channel of the the Wide 
Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST). In the following we 
will focus on the optical data set only, since the subject of the present paper is the accurate 
measure of the stellar PMs, and the larger pixel scale of the WFC3 IR camera detector 
(~ 0.13"pixeR^) makes the secured IR images unsuitable for this purpose. The ACS/WFC 
is made up of two 2048 x 4096 pixel detectors with a pixel scale of ~ 0.05"pixel“^ and 
separated by a gap of about 50 pixels. The total held of view (FoV) is ~ 200" x 200". The 
optical data set consists of 5 x 365 s images in F606W, 5 x 365 s images in F814W, and one 
short exposure per hlter (50 s and 10 s, respectively). The observations were performed on 
August 18th, 2013 and provided an optimal second-epoch data set for PM measures. 


The hrst epoch images (GO-9799, PI: Rich) consist of two deep (340 s) exposures in 
the F606W and F814W hlters, and one short exposure (10 s) in the F814W hlter, also 
acquired with the ACS/WFC. This data set was used to constru c t the deepest optical color- 


magnitude diagram (CMD) of Terzan5 fsee iFerraro et al.l 120091: iLanzoni et ahl 120101 1. The 


images were acquired on September 9th, 2003. Thus, the combination of the two optical 
data sets provides a total time baseline of ~ 9.927 yrs. 


3. RELATIVE PROPER MOTIONS 


The analysis has been performed on _flc images, which have been hat-helded, bias- 
subtracted and corrected for Charge Transfer Efficiency (CTE) losses by the s tandard HST 
calibr ation pipeline (CALACS) adoptin g the pixel-based correction described in lAnderson fc Bedin 
(2010) and lubeda fc Anderson ( 2012 1 


scribed in detail in Massari et al 


The procedure used to derive relative PMs is de- 
fl2013[) . Here we provide only a brief description of the 


main steps of the analysis. The hrst step consists in the photometric reduction of each indi- 
vidual exposure of t he tw o epochs with the public ly available program img2xym_WFC. 09x 10 
Anderson fc King (1200611 : lAnderson et al.l (1200811 . This program uses a hlter-dependent li¬ 
brary of spatially varying PSF models plus a single time-dependent perturbation PSF to 
account for focus changes or spacecraft breathing. The hnal output is a catalog with in¬ 
strumental positions and magnitudes for a sample of sources above a given hux threshold in 
each exposure. Star positions were t hen corrected in each catalog for geometric distortion, 
by means of the solution provided bv lAndersonI (120071) . 
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Figure [T] shows the (mpeoew, ’^Feoew — ’^F 814 w) CMD of Terzan 5 using the second 
epoch data set (Fig. [1]). In particular, for stars fainter than mpeoew — 20.7 mag, which is 
the saturation limit of the deep images (solid grey line in Fig. [T]), all stars detected in at least 



apparent from the figure, the evolutionary sequences of Terzan 5 are clearly distinguishable 
although strongly affected by differential reddening. The main sequence extends for almost 
4 magnitudes below the turn off. A blue sequence is visible at mpeoew < 24.5 mag and 
(?^F 606 W — '^F 8 i 4 w) < 3.6 mag and it remains well separated from the cluster RGB. This 
sequence is likely populated by young disk stars. 


The second step in determining relative PMs is to astrometrically relate each exposure 
to a distortion-free reference frame, which from now on we will refer to as the master frame. 
Since no high-resolution astrometry other than that coming from these data sets is available, 
we defined as master frame the catalog obtained from the combination of all the second- 
epoch single-exposure catalogs corrected for geometric distortion. In this way, the master 
frame is composed only of stars with at least 10 position measurements (5 for each hlter). 
We then applied a counter-clockwise rotation of 91.163° in order to give to the master 
frame the same orientation as the absolute reference f rame defined by the Two Micron All 
Sky Survey (2MASS) catalog fsee lLanzoni et al.ll201Cll) . We then transformed the measured 
position of each star in each exposure into the master frame by means of a six-parameter 
linear transformation. In order to maximizes the accuracy, the transformations have been 
computed using only high signal-to-noise (mpgoew < 24) and unsaturated sources that could 
be considered likely cluster members according to their position in the CMD. Also, in order 
to minimize the effects of GTE- and distrorsion-correction residuals we treated each chip 
separately. At the end of the procedure, for each star we have up to 3 hrst-epoch position 
measurements and up to 12 second-epoch positions on the master frame. To estimate the 
relative PM of each star we computed the median positions in the first and in the second 
epoch by applying a 3cr-chpping algorithm. The difference between the two median positions 
gives the star’s displacement over AT = 9.927 years. Since for stars brighter than mpeoew — 
20.7 one first-epoch and two second epoch positions were available, we adopted the single 
position measured in the hrst epoch and the mean between the two positions measured in 
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Fig. 1.— (mpeoew, wpeoew—™ f814w) CMD of Terzan 5. All the cluster evolutionary sequences 
are broadened because of differential reddening effect. A bright, blue sequence is clearly separated 
from the cluster sequences and it is likely composed of young field stars. The magnitude saturation 
limits of the short and long exposures are marked with the dashed and solid lines, respectively. 


the second epoch. The errors in each direction and within each epoch were computed as: 


where rmsi ,2 is the rms of the positional residuals about the median value, and Ni ^2 is 
the number of measurements. Therefore, the error in each PM-component associated to 
each star is simply the sum in quadrature between hrst- and second-epoch errors: cipj^ = 
a /{<jfy + ((T^)2/AT and = a/T he error associated to the PM of 
the brightest stars measured only in the short exposures were computed by adopting as 
positional uncertainties the typical errors determined in the long exposure catalogs at the 
same instrumental magnitude. 
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Fig. 2.— Uncertainties in the X and Y displacements in units of pixel yr“^. Bright, well measured 
stars have displacement errors typically smaller than 0.007 pixel yr“^ in each coordinate. The 
separated sequence at larger errors which arises for mpeoew < 22.5 is made up of stars measured 
only in the short exposures. 

We repeated the entire procedure several times, in order to reach a stable number of 
stars in the reference list. In performing these iterations, stars were selected on the basis 
of the relative PMs obtained in the previous step. To be conservative, we decided to build 
the final PM catalog taking into account only the 69,425 unsaturated stars measured in at 
least 3 long exposures per epoch and the 899 stars brighter than the deep image saturation 
limit measured in all the short exposures which have a total uncertainty on the displacement 
smaller than 0.2 pixels. All the sources which saturate in the latter images have been 
excluded. The typical error as a function of magnitude is shown in Figure [2l For well- 
exposed stars it is smaller than 0.007 pixel yr“^ in each coordinate. Faint stars or stars with 
only few measurements show larger errors, still typically smaller than 0.02 pixel yr“^ (see 
the separated sequence in Fig. |2]). 














7 


3.1. The Vector Point Diagram 


We converted the PMs into nnits of masyr“^by mnltiplying the measnred displacements 
by the pixel scale of the master frame (0.05"/pixel) . Since the master frame has been 
oriented according to the equatorial coordinate system, the X PM-component corresponds 
to that projected along (negative) Right Ascension (—/iQ,cos5), while the Y PM-component 
corresponds to that along Declination The output of this analysis is summarized in 

Figure [3l where we show the vector point diagram (VPD) for all the stars of the final PM 
catalog. The first clear feature is that more than 70% of the stars are distributed within the 
innermost 1.5 masyr“^ from the origin of the VPD (see histograms in Fig. [3]). The remaining 
fraction of stars describes a sparser and asymmetric distribution out to about 10 masyr“F 
In order to highlight these features, we selected only stars with mpeoew < 24, which typically 
have the most accurate PMs, and plotted them as red points in Figure |3l Their distribution in 
the VPD clearly shows at least two components. One is a symmetric distribution centered 
around the origin, with PM weighted mean values HaCosS = 0.02 ± 0.02 rnasyr”^ and 
Us = —0.01 ± 0.02 masyr“^. The other is an asymmetric structure approximately centered 
around the coordinate {^aCos6 ~ —1.8, ~ —3.5) masyr”^ in the VPD. The location of 

these two components in the CMD clearly reveal their nature (see Figure 0]). In fact, the stars 
of the first component (shown as blue dots in the VPD) describe the cluster evolutionary 
sequences in the CMD (central lower panel), with a small degree of contamination left. On 
the other hand, the stars belonging to the asymmetric component (red dots) correspond to 
the blue plume in the CMD (right-lower panel), which is essentially populated by young 
disk stars in the foreground of Terzan 5. This is further c onhrmed by the comparison with 
the prediction of the Besangon model flRobin et al.l l2003l) with only young {tags < 7 Gyr) 
Galactic disk stars for a held centered at the coordinates of Terzan 5 and covering the same 
FoV as that of the AGS/WFG (see Figure [5]). 


3.2. Check for systematic errors 


Many systematic uncertainties may affect the measure of relative PMs (see iBellini et ah 


2014J for a detailed description). It is therefore important to check that the measured PMs 
of our dual catalog do not suffer from such systematic effects. Since we have only two 
epochs of observations and a relatively small number of images per epoch (only two long 
exposures were acquired in the hrst epoch) the full control of all the possible systematics is 
not achievable. Nonetheless, we have carefully checked for the presence of any systematics 
that our data sets allow us to verify. 


First of all we looked for any possible chromatic-induced systematics by checking that 








Fig. 3.— Vector point diagram (VPD) of the PMs measured for 70,324 stars (black dots) in the 
direction of Terzan 5. Their distribution in the Right Ascension and Declination PM-component 
axes is shown in the histograms in the bottom and right panel, respectively. PMs measured for 
unsaturated stars with mpeoew < 24 mag are shown as red dots. At least two components are 
visible: the first showing a symmetric distribution centered around the origin and an asymmetric 
structure roughly centered at (-1.8, -3.5) masyr“^. 


no trend between PMs and color exists for our sample. Figure [6] shows the results of this 
test. By selecting stars in a wide magnitude range (23 <mpeoew < 26), we computed the 3-a 
clipped average value of our PMs (for both the spatial components) in color bins of 0.04 mag 
(red hlled circles). The errorbars are within the size of the symbols. Clearly, these values 
describe a flat relation (shown as a red dashed line). The best least square linear £t to these 
data gives a null angular coefficient within the uncertainty for both the PM components. 
Thus, we can safely conclude that no trend with color exists for our measurements. 

We repeated a similar test looking for possible trends with the observed magnitude in 
the range 19 < mpgoew < 26, thus to exclude extremely faint sources that have intrinsically 
larger PM errors. Figure [7] shows the result of this test. Also in this case, the best linear 
£t to the 3-cr clipped average PM values computed in bins of 0.3 mag is compatible (within 
the errors) with a null angular coefficient line, thus demonstrating that no significant trend 
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Fig. 4.— Upper panels: VPD of all stars brighter than m.Fme,w = 24 mag. The PMs of stars 
belonging to the symmetric component of likely members are highlighted in blue in the central 
panel, those of stars belonging to the asymmetric component centered at (—1.8, —3.5) masyr“^ 
are plotted in red in the right-hand panel. Lower panels: CMDs described by the PM-selected 
stars in the corresponding VPDs. When likely Terzan 5 members are selected, only the cluster 
evolutionary sequences are visible, with the exception of few residual contaminating stars. On the 
other hand, the population belonging to the asymmetric component appears to be dominated by 
young foreground disk stars along the blue plume. 


is found also with magnitude. 

Finally, we also carefully checked the existence of spurious trends with the loc a tion o f 
stars on the detector. To do so, we followed the method described in iBellini et al.l f)2014j) . 
by plotting maps of the measured PMs as a function of the position on the sky, where to 
each star we associated the average motion of the closest 100 neighbors. Overall we found a 
homogeneous distribution in both the PM components. 


The tests discussed above demonstrate that the derived PMs do not suffer from signif¬ 
icant systematic effects. Also, the use of available GTE corrected images for both epochs 
(see SectlH]) should have minimized the effect of GTE losses on the PM estimate. While we 
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Fig. 5.— Comparison between the observed ACS CMD and that predicted by a simulation of the 
Besangon model including only Galactic disk stars younger than 7 Gyr. Such a comparison clearly 
demonstrates that the bluer sequence in the CMD of Terzan 5, already identified as composed of 
field stars by the measured PMs, correspond to the MS of foreground disk stars. 

cannot exclnde that other systematics can affect ont measnrements, we can firmly conclude 
that these are the best PMs we can measure with the available data set. 


4. RESULTS 

We used the measured relative PMs to “clean” the optical (mpeoeW) ~ ''^F 814 w) 

CMD of Terzan 5 described in Section [2J We dehned as likely member stars all the sources 
located within a distance smaller than 1.5 masyr“^ in the VPD plane. As shown in Figure El 
this appears to be a reasonable assumption, since the bulk of the stars lie within this limit. 
The CMD obtained from such a selection is shown in Figure El Member stars are plotted 
as black dots, while non members stars are shown in red. The selection applied leaves in 
the CMD only stars clearly belonging to the cluster evolutionary sequences, while excluding 
most of the outliers. A small degree of contamination is still present because the distribution 
of held stars (mainly bulge stars) in the VPD overlaps that of Terzan 5 members. However, 
we can conclude that the PM analysis performed is efficient in decontaminating the CMD 
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Fig. 6.— Upper panel: Y-component of the relative PMs vs observed colors. Red filled circles 
represent the 3-a clipped mean PMs measured in color bins of 0.04 magnitudes. Stars rejected by 
the clipping algorithm are shown in grey. All the mean values lie on the no-correlation line (red 
dashed line), thus ecluding any chromatic-induced systematics. Lower panel: same as above, for 
the PM X-component. 


from foreground and background sources. We ca n therefore us e this se lection to ass e ss the 


mern b ership of the sub-popu lations discussed bv iFerraro et al.l 020091) : lOriglia et al.l (12011 


2 OI 3 I ): iMassari et al.l f 2014bl ). 


4.1. MAD Infrared CMD 


The two major su b-populations in Terzan 5 were identified in the {K, J — K) CMD 
bv iFerraro et al.l (120091) . Their radial distributions are incompatible with that expected for 
background held stars, and the measured line of sight velocities are consistent with the 
systemic radial velocity of Terzan 5. However, the membership of the two populations has 
been questioned (see, e.g., IWillman fc Straderll2012l) . To address this concern and to assess 
their cluster membership from proper motions, we cross-correlated the MAD catalog with 




















12 



Fig. 7.— Upper panel: same plot as in Figure[6l but in terms of observed magnitude. Mean values 
have been computed in bins 0.3 mag wide. Also in this case, no trend is found between PMs and 
magnitude. Lower panel: same as above, for the X-component of PMs. 


the PM catalog. The result is summarized in Figure [H The left panel shows the IR CMD 
of Terzan 5 after the PM-based membership selection (black dots). The sources recognized 
as contaminating held stars based on the value of their PMs are also plotted in red. The 
decontaminat e d CM D clearly exhibits the two red clumps (RCs) originally discovered by 


Ferraro et al.l ( 20091) . 


Figure [TO] highlights the CMD and VPDs for sub-samples of stars properly selected 
in the two RCs. As can be seen, the two PM distributions appear quite symmetric, both 
showing a small (0.5 rnasyr”^) dispersion around the origin. While these PMs cannot be 
used to reveal possible intrinsic kinematical differences between the two populations, we can 
solidly conclude that these stars are all members of Terzan 5. 


Overall the number of contaminants in the IR CMD (selected as stars outside a distance 
of 1.5 mas yr“^ from the center of the VPD) is much smaller than what observed in the optical 
plane. This is because the MAD photometry corresponds to a smaller FoV (with respect 
to the ACS one) and to the very central region of the system, where the cluster population 
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Fig. 8.— Left panel: optical CMD of Terzan 5, with likely member stars shown as black dots, and 
field sources (excluded by the PM-based selection) plotted in red. The measured relative PMs are 
efficient in decontaminating the CMD and only cluster evolutionary sequences survive the selection 
criterion. Right panels: magnitude-binned VPDs for all the stars in the optical catalog. Each bin 
has a size of 2 mag. Sources are color coded as in the left panel. 

is expected to dominate. It is also worth noticing that both the RC sub-samples are much 
more centrally concentrated than the likely Galactic held contaminants selected at the same 
magnitude level (see the red, blue and grey lines, respectively, in FiglTOl). According to a 
Kolmogorov-Smirnov test performed on the three samples, the probability of PM-rejected 
stars to have been extracted from one of the two RC populations is always smaller than 
10“^. This further conhrms that stars rejected according to our PM-membership selection 
are held objects, homogeneously distributed in the sampled FoV. 

In conclusion, the relative PMs measured in this work and applied to the MAD CMD def¬ 
initely demonstrate that the two major sub-populations photometrically discovered in Terzan 
5 are both members of the system. 
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Fig. 9.— Left panel: IR CMD of Terzan 5 obtained from MAD observations. The PM-selected 
member stars are shown as black dots, while the sources having discordant PMs are marked in red. 
The presence of two distinct red clumps is well evident even after the PM-based selection. Right 
panels: magnitude-binned VPDs (each magnitude bin has a size of 3 mag). Stars are color-coded 
as in the left panel. 


4.2. Spectroscopic targets 


In iMassari et al.l fl2014bl) we presented the distribution of the iron abundance for the 
largest sample of stars in Terzan 5. In that work, target membership was inferred from radial 
velocities and the fraction of contaminating stars was estimated statistically fr om the chemi¬ 


cal an d kinematical characterization of the held surrounding the system itself (IMassari et ah 


2014a|). According to that work, in the innermost 200"the expected fraction of contaminating 


held stars is about 16% (18 out of 114 stars). 


The measure of relative PMs gives the great opportunity to verify the membership of 
each single star in the ACS FoV (instead of adopting a statistical approach only) and to 
check whether the estimate of the contribution from contaminating stars is reliable. To this 
aim, we cross-correlated them with the iron abundance catalog. We found only 42 out of 
135 targets in common, because the majority of the spectroscopic targets are either bright 












15 



(•^ ^)mad 


r" 2 

I 

>% 

« 0 

B 

^ -2 

^ 2 

u 

>^ 

a 0 

B 

^ -2 


/i„cos((5) [mas yr“‘] 

1 

0.8 

2 0.6 
z 

Z 0.4 
0.2 
0 

0 10 20 30 40 

dist[arcsec] 




2 0-2 


Fig. 10.— Left panel: IR MAD CMD of Terzan 5 zoomed in the RC region. Faint-RC stars 
are plotted in blue, while bright-RC stars are plotted in red. Stars rejected according to their 
PMs are shown in grey. Upper- and central-right panels: VPDs of the bright-RC and faint-RC 
sub-samples, respectively: the stars of both populations lie within the 1.5 masyr“^ circle adopted 
as membership selection criterion. No clearcut difference between the two distributions is visible. 
Lower-right panel: cumulative radial distributions of bright-RC, faint-RC and likely non-members 
stars selected in the IR CMD (red, blue and grey lines, respectively). Clearly, non-member stars 
are less concentrated than cluster members, as is expected for field sources. 


RGB stars, which saturate even in the short exposures, or stars located outside the ACS 
FoV. The location of these targets in the VPD is shown in Figure [HI where blue circles 
mark stars with —0.6 <[Fe/H]< 0, while red circle correspond to more metal rich objects 
(0 <[Fe/H]< 0.5). Clearly, with only a few exceptions (4 targets), all spectroscopic targets 
appear to be members of the system, having PMs (and radial velocities) well within the 
characteristic cluster distribution. Also four stars belong ing to the newly discovered metal 


poor c omponent of Terzan 5 (with [Fe/H]~ —0.8, see lOriglia et ah! 120131: iMassari et ah 


2014bh have measured PMs and are shown in the hgure (black triangles). As apparent, they 
are also members of Terzan 5. Finally, the t wo magenta hll e d squa res correspond to the 
super metal-rich stars ([Fe/H]> 0.5) found in Massari et ah f 2014b[) and suggested to be 
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non-members. Indeed, their PMs are either outside the distribution of member stars, or just 
at its edge. 


Overall, 6 out of 42 spectroscopic targets appear to be held stars. Such a fraction 
corresponds to the 14% of the sample, in very good agreement with the statistical estimate 
obtained from our previous spectrosco pic studies. Therefore , these results conhrm that the 
iron distribution of Terzan 5 shown in iMassari et al.l fj2ni4m is made up of genuine cluster 
members, showing a huge internal iron spread of more than 1 dex. 



Fig. 11.— VPD for the spectroscopic targets discussed in IMassari et al.l (|2014bl ) in common with 
our PM catalog. Blue circles are stars with —0.6 <[Fe/H]< 0, while red circles have 0 <[Fe/H]< 0.5. 
These represent to the two main populations of Terzan 5. Black filled triangles correspond to stars 
belonging to the newly discovered third, metal-poor population. All the three components are 
clearly genuine members of the system, being located well within the PM distribution of cluster 
members. The two magenta filled squares represents the super metal-rich stars of the spectroscopic 
sample and, at odds with the other sources, they seem not to belong to Terzan 5. 
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4.3. Absolute proper motion 

One of the most intriguing application of PM measures is the determination of the 
absolute motion of a stellar system. The best way is to relate the mean motion of all the 
stars in the systemic to the absolute reference frame dehned by very distant o bjects, such as 


background galaxies or quasars, which appear as zero-motion sources (see, e.g., Dinescu 

et al. 

1999; 

Mahmu 

d & Anderson 

2008: Bellini et al. 

2010 ; 

Lenine et al. 

2011 ; 

Sohn et al. 

2 OI 2 I; 

Massari et al. 

2013 

). Unfortunately, because of the high stellar density and the very large 


extinction in the direction of Terzan 5, we have not been able to detect such objects in the 
ACS FoV. 

We then cro ss-correlated our PM catalog with many public catalogs of stellar PMs, 
such as NOMAD f Zacharias et ah 2005 ). the Guide Star Catalog (GSC) version 2.3, UCAC4 


fIZac 


larias 


1998 


et ah ( 20121) j a nd the Yale/San Juan Southern Proper Motion catalog f Platais et al 


Dinescu et ah 


Il997h 


However, since these catalogs have PMs measured only for very 
bright sources, we found just fours stars in common with our PM sample, all having PM 
uncertainties larger than 8 masyr“^ in the public data sets. Hence, no meaningful result 
about the absolute PM of Terzan 5 can be obtained from this kind of analysis. 


' Ortolan! et al. 

2011 

, see also 

Rossi et al. 


20151) . where the authors anchor the motion of the GC HPl to that of the underlying 
bulge population. The bulge PM is the composition of its internal kinematics and the re¬ 
flex motion of the Local Standard of Rest (LSR). In the direction o f Terzan 5, the radia l 
velocity distribution of bulge stars peaks at Urad = 21.0 ± 4.6 km s“^ flMassari et al.ll2014al) . 


Assuming that bu lge stars are at the same distance of t he Galactic ce nter (here we adopt 


Ghez et al. 

2008; 

Ortolan! et al. 

2011 

, but see also 

van der Marel et al. 


2 OI 2 I for a complete overview of the topic), the tangential component of the bulge velocity is 
Vta.n = 'yradsin(/)/cos(/) = 1.41 ±0.31 km s“\ or in Gal actic coordinat e s = 0.05 ±0.01 


mas y r 


-1 


Since t h e bulge shows cylindr i cal ro tation [Howard et al.l (120081); iKunder et al. 


(2012); Ness et al. f 2013b[ l: Zoccali et al. (2014), we can assume = 0 masyr“L By 

summing these values to the motion of the LSR (ui cosjb) = 6.10 ± 0.25 mas y r~^, = 0 

masyr”^, corresponding to ulsr = 243 ± 10 km s“^; see lOrtolani et al.l 1201 ll) . we obtain 


the total motion of the bulge at Terzan 5 coordinates: (/x; cos(6), /i;,)buige = (6.15 ± 0.25, 0) 
mas yr“^. 


As a second step for determining the PM of Terzan 5 relative to that of the bulge. 


^Given the large number of measured PMS, this kind of analysis usually does not require the accuracy 
needed to study the internal kinematic of a system. 
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we excluded foreground disk contaminants belonging to the blue plume in the CMD by 
selectining a sub-sample of stars with (mFeoew-^ipsuw) > 3.9 and mpeoew < 24.5, which do 
not saturate even in the long exposures. We also applied a selection in terms of PM errors, 
by excluding stars with uncertainties larger than 0.3 rnasyr”^ in the two PM components. 
We then defined as Terzan 5 members all the stars located within 2masyr“^ from the origin 
of the VPD. 

To select the sample of bulge stars w e foll owed an iterative procedure similar to that 
described by [Anderson fc van der Marell (12010!) for the determination of the center of u 
Centauri. After the exclusion of all Terzan 5 members, we computed the hrst guess mean 
motion of bulge stars. We then excluded all objects contained within a circle of 2masyr“^ 
radius positioned at the same distance as the “Terzan 5 circle” from the hrst guess mean, 
on the symmetrically opposite side in the VPD. A new value of the mean motion has been 
evaluated from this sample, and the procedure has been repeated until convergence was 
reached. Once the Terzan 5 and the bulge samples are dehned, we determined their weighted 
mean motion with a 3-a clipping procedure that, after the last step, included 3853 sources 
for the Terzan 5 sample and 797 for that of the bulge (black dots in Figure IT^. We computed 
the uncertainty in the weighted mean motions as the standard error-in-the-mean, i.e., the 
dispersion of the surviving stars around the mean PM, divided by the square root of their 
number. 

The resulting motion of Terzan 5 relative to the bulge (see the blue arrow in Figure IT^ 
expressed in Galactic coordinates (see lOrtolani et ahll201ll for the details of the transforma¬ 
tions) turned out to be (/i/cos(6), /i 6 )T 5 -buige = (0.26 ± 0.10, 0.83 ± 0.12) masyr”^ , where 
the final uncertainties are the sum in quadrature of those computed for the two separated 
samples. By subtracting the motion of the bulge quoted above, we finally obtained the ab¬ 
solute PM of Terzan 5: (/i;cos(6), ^J,b)T 5 = (—5.89 ± 0.27, 0.83 ± 0.12) masyr“^, where the 
total uncertainties are the sum in quadrature between that coming from the motion of the 
bulge and that associated to the motion of the cluster relative to it. This is the hrst absolute 
PM estimate ever obtained for Terzan 5. 


4.4. The cluster orbit 


On the basis of the absolute PM estimate derived in the previous section, we performed 
a numerical integration of the orbit of Terzan 5 in the Gal actic potential. We used the 3- 


component (bulge, disk and halo) axisymmetric model from lAllen &: SantillanI fll99ll) . which 
has been extensively used and discussed in the liter 
vironmental effects on Galactic stellar systems (e.g 


Allen et al.l 

2006; 

Montuori et ah 

2007 
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Fig. 12.— VPD of the sample of stars used to estimate the absolute PM of Terzan 5. The black 
points correspond to the bulge stars surviving the last step of the selection procedure and effectively 
used as reference; the gray dots mark the Terzan 5 member stars. The mean motion of the bulge 
sample is marked with a red cross. The PM of Terzan 5 relative to the mean bulge motion is 
marked with a blue arrow. 


Ortolani et al.l l201ll: iMoreno et al.l 120141: IZonoozi et al.l l2014i) . thanks to its relative sim¬ 


plicity and fnlly analytic natnre. We adjnsted the various model parameters to make the 
rotation velocity curve match the value of 243 km s“^ measured at the Solar Galactocentric 
distance of 8.4 kpc (see above). The cluster PM were reported in the cartesian Galactocen¬ 
tric reference frame, resulting in a velocity vector {vx,Vy,Vz) = (—60.4 ± 1.3, 85.7 ± 11.1, 
35.0 ± 6.0) km s“^ at the position {x,y,z) = (—2.51 ± 0.30, 0.39 ± 0.02, 0.17 ± 0.01) kpc. 
We adopted the convention in which the X axis points opposite to the Sun (i.e., the Sun 
position is (—8.4, 0,0)). The orbit was then time-integrated backwards for 12 Gyr, starting 
from the given initial (curr ent) conditions and using a 2nd order leap-frog integrator (e.g. 
Hockney fc Eastwood 19881) with a rather small and constant timestep (~ 100 kyr, corre¬ 


sponding to ~ 1/300 of the dynamical time at the Sun distance, computed as (i?g/GMg)^/^, 
where Mg ~ 1.7 x IO^^Mq is a characteristic Galactic mass parameter). During the > 10^ 
timesteps used to describe the entire orbit evolution, the errors on the conservation of both 
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the energy and the Z component of angnlar momentnm were kept under control and never 
exceeded 1 part over 10^ and 10^^, respectively. To take into account the uncertainties on 
the kinematic data, we generated a set of 1000 orbits starting from phase-space initial condi¬ 
tions normally distributed within a 3cr range around the cluster velocity vector components 
and current position, with a being equal to the quoted uncertainties on these parameters. 
For all these orbits we repeated the backward time-integration. The probability densities of 
the resulting orbits projected on the equatorial and meridional Galactic planes are shown 
in Figures (TT] and dU respectively. Darker colors correspond to more probable regions of 
the space, i.e., to Galactic coordinates crossed more frequently by the simulated orbits. As 
apparent, the larger distances reached by the system during its evolution are R = 3.5 kpc 
and \Z\ = 1.6 kpc at a 3cr level of signihcance, which roughly correspond to a region having 
the current size of the bulge. 


In a forthcoming study, we plan to quantify the possible influence of a central non ax- 
isymmetric component (the bar and/or the spiral arms). Indeed, although triaxial potentials 
are expected to facilitate the onset of very radial orbits (because of the lack of conservation 
of any angular momentum component) with a p ossible prolon g ation of the system orbital 
motion in much outer regions, a recent study by iMoreno et al.l (120141) demonstrated (using 
a sophisticated and realistic model) how this intuitive picture could be wrong, due to the 
complex interplay between the gravitational potentials generated by the various components. 
The relatively small perigalactic distance (< 500 pc) that Terzan 5 could have achieved dur¬ 
ing its orbit could imply a significant mass loss due to tidal erosion and thus represents a 
further strong motivation for studying in more detail the dynamical history of this system. 


Based on the current results, we can conclude that the integration of Terzan 5 orbit 
performed by means of this simple axisymmetric Galactic static model suggests an in-situ 
formation for the cluster wi thin the Galaxy, ra ther than an external accreted origin. This 
supports the interpretation (iFerraro et al.l 120091 1 of this system as the remna nt of one of the 


massive stellar clum ps that may have contributed to form the Galactic bulge fjNoguchilll999 


Immeli et al.l 120041) . 


5. CONCLUSIONS 

By using high-resolution AGS/WFG images, we have measured relative PMs in the 
direction of the stellar system Terzan 5. We have been able to separate cluster members 
from foreground disk stars and background bulge contaminants. The PM-selected GMD 
clearly conhrmed that the two main populations with distinct iron content discovered in 
Terzan 5 are genuine members of the system, thus putting to rest this issue which had 
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Fig. 13.— Probability density in the eqnatorial Galactic plane of 1000 simnlated orbits of 
Terzan 5 time-integrated backwards for 12 Gyr. Darker colors correspond to larger proba¬ 
bilities. 



Fig. 14.— As in Fig. [131 bnt for the orbit projection in the meridional Galactic plane. 
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remained a matter of debate. Moreover, the measured PMs demonstrate that also the 
third, metal-poorer and a-enhanced population clearly belongs to the system. Our findings 


therefore conhrms t hat Terzan 5 is, together with u Centauri flNorris &: Da Costa 


Pancino et al.ll2000l: lOriglia et al.ll2003l: I.Tohn.son fc Pilachowskill2010l: IVillanova et al 


1995 


291411 . 


the stellar system with the largest internal iron spread in the Galaxy. 


Finally, by comparing the motion of Terzan 5 members to that of bulge stars, we 
estimated for the first time the absolute PM of this system, finding (/i;cos(&), /i 6 )T 5 = 
(—5.89 ±0.10, 0.83 ±0.12) masyr”^. The backward integration of its orbit, computed under 
the assumption of an axisymmetric Galactic model, provides another indication that sup¬ 
ports the scenario according to which Terzan 5 is the remnant of one of the primordial stellar 
clumps that may have contributed to form the galactic Bulge, and not an object accreted 
from outside the Milky Way. 


We thank the anonymous Referee for his/her useful comments which helped us to im¬ 
prove our paper. The authors are also grateful to Dr. Lucie Jilkova for her valuable sugges¬ 
tions on the cluster orbit integration. This research is part of the project Cosmic-Lab (web 
site: http://www.cosmic-lab.eu) funded by the European Research Gouncil (under contract 
ERG-2010-AdG-267675). D.M. thanks the Kapteyn Astronomical Institute for hospitality. 
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The techniques applied i n the present work have been developed in the context of the HST- 
PROMO collaboration!^ (j van der Marel et ah! I 2 OI 4 I: iBellini et ah! I 2 OI 4 I) , which aims at im¬ 
proving our understanding of the dynamical evolution of stars, stellar clusters and galaxies 
in the nearby Universe through the measurement and interpretation of PMs. 
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